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Abstract. Eclipsing binary stars are gratifying objects because of their unique 
geometrical properties upon which all important physical parameters such as masses, 
radii, temperatures, luminosities and distance may be obtained in absolute scale. 
This poses strict demand on the model to be free of systematic effects that would 
influence the results later used for calibrations, catalogs and evolution theory. We 
present an objective scheme of obtaining individual temperatures of both binary 
i-C ' system components by means of color-index constraining, with the only requirement 

^ , that the observational data-set is acquired in a standard photometric system. We 

show that for a modest case of two similar main-sequence components the erroneous 
^\ \ approach of assuming the temperature of the primary star from the color index 

yields temperatures which are systematically wrong by ~ 100 K. 
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^ ! 1. Introduction 

o 

\ Eclipsing binary stars (EBs) are renowned for their typical geometri- 

cal layout and well-understood underlying physics. Solving the inverse 
^ . problem for detached EBs has become an indispensable tool for obtain- 

^ \ ing absolute dimensions of individual stars (their masses, radii, tem- 

^ ■ peratures, and luminosities). The impact is very broad: EBs are used 

"to . to establish color calibrations (Harmanec, 1988; Flower, 1996; Pop- 

per, 1998), evolutionary properties based on exact coevality of both 
components (Ribas et al., 2000; Lastennet and Valls-Gabaud, 2002) 
^ \ and relations between M-L-R-T quantities for different spectral types 

• and luminosity classes (Gorda and Svechnikov, 1998). It is thus very 

important to model EBs as objectively as possible, otherwise it is likely 
that systematic errors may creep in and mislead any theoretical models 
built on those results. In this paper we concentrate on a particular 
issue of objectively determining effective temperatures of individual 
components in an EB system. 

The common practice of determining effective temperatures of indi- 
vidual components is to assume the temperature of one star and obtain 
the temperature of the other star by the model. The assumed tempera- 
ture is usually obtained from observed color-indices at quarter-phase by 
using a color calibration or by spectral analysis and comparison against 
synthetic spectra database. While this approach may be adequate for 
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Table I. Physical parameters of the F8 V-Gl V test binary star. Linear (x) and 
non-linear (t/) coefficients of the logarithmic limb darkening law for Johnson B and 
V passbands, taken from van Hamme (1993). 



V 



¥ 



Parameter [units] 



Po [days] 
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Vj [km s ] 

Teff [K] 

L [Re] 

M [Rq] 

R [Re] 

n [~]'^ 

log (g/gu) 

Xb 

Vb 

XV 

yv 



F8 V 



Binary 



6200 
2.100 
1.239 
1.260 
5.244 
4.33 
0.818 
0.203 
0.730 
0.264 



1.000 
5.524 
0.831 
85.000 
15.000 



Gl V 



5860 
1.100 
1.030 
1.020 
5.599 
4.43 
0.833 
0.158 
0.753 
0.242 



x/a 



°'fi is potential as defined in Wilson (1979). 
'go = Icms^^ makes g/go dimensionless. 



components having significantly different luminosities, (so tfiat the lu- 
minosity of one star is equal to the system luminosity for all practical 
purposes), it is very inadequate for components with similar luminosi- 
ties and introduces systematical effects into model solutions. When 
dealing with photometric observations that are acquired in a standard 
photometric system (i.e. Johnson, Cousins, Stromgren, . . . ), effective 
temperatures of both components may be disentangled from two or 
more photometric light curves objectively, without any assumptions. 
We build this argument and demonstrate its implications on modeling 
in the following Sections. 

2. Simulation 

To evaluate the uncertainty of determining a primary star temper- 
ature a priori, we built a synthetic binary star model. Testing the 
method against a synthetic model may seem artificial, but the obvious 
advantage of knowing the right solution is the only true way of both 
qualitative and quantitative assessment. 

Our synthetic binary consists of two main-sequence F8 V-Gl V 
components with their most important orbital and physical parameters 
listed in Table I. It is a partially eclipsing detached binary with only 
slight shape distortion of both components (i?i, pole /-Ri, point = 0.974, 
-R2,poie/-R2,point = 0.979). Light curves are generated with PHOEBE (Prsa 
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and Zwitter, 2005) for Johnson B and V passbands in 300 phase points 
with Poissonian scatter cjlc = 0.015 and quarter-phase magnitude 
my = 10.0. Passband transmission curves were taken from ADPS 
(Moro and Munari, 2000). 

The simulation logic is as follows: we take Kurucz's spectral energy 
distribution (SED) function for both components from the spectra 
database compiled by Zwitter et al. (2004), doppler-shift them to the 
reference phase, sum their SEDs weighted by their corresponding lu- 
minosity and convolve the sum with the passband response functions 
of Johnson B and V filters. This enables us to compute the B-V color- 
index, which needs to be zero-corrected with respect to spectral type 
AO (adopted temperature 9560 K). Once this is done, we have obtained 
a synthetic prediction of the color index. Alternatively, the same result 
may be obtained by using color calibrations, e.g. Flower (1996) 's cor- 
rected tables given by Prsa and Zwitter (2005). Color indices may now 
be readily transformed to the ratio of passband fluxes J^^/J^f. 

Once the light curves are built so that the color index is preserved, 
we follow the color-index constraining principle described by Prsa and 
Zwitter (2005) to assess the implications of erroneous assumed value of 
Teffi- 3. Results 



Our analysis brings us to the following conclusions: 

1. Individual temperatures without the imposed color-index constra- 
int are fully correlated, as is well known. The reason for this is the 
degeneracy between the temperatures and passband luminosities 
L\: local effects of slightly different temperatures on the shape of 
the light curve are very small, hardly noticeable (they scale roughly 
with the ratio of both temperatures). However, even slightly differ- 
ent temperatures cause an additional effect which is an order of 
magnitude larger: they scale both light curves, thus changing the 
color index. If this change is not constrained, the model relies only 
on secular effects, which are almost always buried in data noise and 
parameter degeneracy. Without color-constraining there is thus no 
practical way to pinpoint a location in the T^si-Tetei cross-section 
that corresponds to the right solution. 

2. Unless the luminosity of one component is overwhelming, the as- 
sumed value of the temperature is doubtful. Fig. 1 shows a range of 
assumed values for Teffi and the color indices yielded by the model 
that is not color-constrained. Since shapes of the light curves deter- 
mine the ratio of luminosities of individual stars, it has to remain 
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Figure 1. Color index dependence on the assumed value of the primary star temper- 
ature. If color-constraining is not applied, all depicted points may be satisfied by the 
model. If we are to assume the value of Tcfri from the color index at quarter-phase, 
it would yield Tefii = 6080 K, a value which is ~120K off from the correct value. 
If on the other hand we apply color-constraining, we immediately get a handle on 
which color index corresponds to the temperatures yielded by the model. Cross-hairs 
denote the solution obtained by color-constraining, which gives Tefii = 6187 ± 32 K. 



approximately constant. To achieve this with changing tempera- 
tures, the model adapts passband luminosities (first order effect) 
and stellar radii (second order effect) accordingly. The problem 
now becomes immediately evident: if we are to assume a primary 
star temperature from the quarter-phase color index, it would be 
roughly T^a|sumed ^ 6080 K, around 120 K off just because the lu- 
minosity of the secondary star is not negligible. In principle, the 
solution could then be manually re-iterated by computing the pre- 
dicted color-index from the solution, comparing it to the observed 
color index and then modifying the assumed temperature to get 
closer to the right solution. This is tedious and prone to subjective 
considerations, it should thus be avoided. 



3. Adopting the color-constraint means keeping the ratio of passband 
luminosities L\ /Lf accordant with the observed color index. This 
means that only LY is adjusted to reproduce the V light curve and 
Lf (and other L\ for the remaining passbands) is then computed 
from this constraint. The degeneracy between the temperatures 
and the passband luminosities is now broken: we obtain a straight 
horizontal line over Fig. 1 which best represents the color-index. 
Thus, we obtain the primary temperature without any assumptions, 
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Figure 2. Convergence tracer for the Tcti'2/Tcsi-^2 cross-section. Cross-hairs denote 
the correct solution. Panels (a) to (d) correspond to the solution obtained by Nelder 
& Mead's downhill Simplex (NMS) for 0, 1, 2 and 3 paxameter kicks. See Prsa & 
Zwitter (2005) for details on interpreting convergence tracers. 



and the secondary temperature follows directly from the model. 
Applying this to our test binary, we obtain individual temperatures 
Teffi = 6 187 K ± 32 K and Teff2 = 5 868 K ± 31 K follows from the 
model. Comparing these values to true values T\ = 6200K and 
T2 = 5 860 K is very encouraging. 

4. Finally, let us consider the correlation between the temperatures 
and the radii of the binary system components. We have men- 
tioned earlier that this is a second order effect, but it is still very 
important to quantify the correlation between those two param- 
eters. Since the radii come into the model via surface potentials 
ri (defined in Wilson, 1979), we need to consider the correlation 
between the ratio of temperatures T2IT1 and surface potentials. 
Fig. 2 shows the convergence tracer for the T2/T1-Q2 cross-section: 
tracing parameter values from each starting point, iteration after 
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iteration, all the way to the converged solution. Attractors - regions 
that attract most convergence traces - within these cross-sections 
reveal parameter correlations and degeneracy. Panels denoted with 
(a) to (d) represent three subsequent parameter kicks, which are 
introduced and thoroughly explained by Prsa and Zwitter (2005). 

4. Discussion 

We have demonstrated how severe the implications of assuming a 
primary star temperature from the color index may be if extra care is 
not taken to submit the acquired solution to manual re-iteration. The 
effect is already very pronounced for a modest case of two similar main- 
sequence components and it would be significantly augmented e.g. for 
a hot main-sequence-cool giant system, where both stars have similar 
luminosities; in those cases the published solution may be off by as much 
as ~1000K. By having the color-constraining method in PHOEBE, we 
hope to raise awareness of this issue in future solution-seeking analyses. 
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